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We discuss the general operator product expansion of a non-forward unequal mass virtual 
compton scattering scattering amplitude. We find that the expansion now should be done 
in double moments with new moment variables. There are in the expansion new sets of 
leading twist operators which have overall derivatives, and they mix under renormalization. 
We compute the evolution kernels from which the anomalous dimensions for these operators 
can be extracted. We also obtain the lowest order Wilson Coefficients. In the high energy 
limit we find the explicit form of the dominant contributing anomalous dimensions and solve 

■ the resulting renormalization group equation. We find the same high energy behavior as 
I/") ' indicated by the conventional double leading logarithmic analysis. 
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ON 

O I. INTRODUCTION 

qj i Recently there is much interest in the analysis of non-forward scattering processes, for example, Deeply 

Virtual Compton Scattering (DVCS) [QJ|] and hard diffractive electroproduction of vector mesons H ^| in 
Deeply Inelastic Scattering (DIS). 

Ji (l]J^] proposed that one can obtain from DVCS information on the Off-Forward Parton Distributions 
(OFPD) which in this case contain new information on long distance physics. He studied their evolution and 
sum rules (in terms of form factors) and obtained certain estimates at low energy. Radyushkin also studied 
the scaling limit of DVCS and generalized the discussion to hard exclusive electroproduction processes 
§|. The non-perturbative information is incorporated in his double distributions F(x, y; t) and non-forward 
distribution functions F^(X;t). He discussed their spectral properties, the evolution equations they satisfy, 
their basic uses and general aspects of factorization for hard exclusive processes. A third parametrization 
for the non-perturbative information was proposed by Collins et al. J^|^] . They performed a numerical study 
of their non-diagonal parton distributions in leading logarithmic approximation, in which they found the 
nondiagonal gluon distribution X2G(xi,X2,Q 2 ) can be well approximated at small x by the conventional 
gluon density xG(x, Q 2 ). 

In this paper, we formulate a general operator product expansion description for a generic non-forward 
unequal mass virtue compton scattering amplitude(see Fig. [l]). Because of the non zero momentum transfer 
?' = p — p' from the initial state proton to the final state proton, the scaling behavior must be different 
from that in a forward case because of the new kinematic degree of freedom (the non- forwardness). Indeed 
we will show that compared with forward scattering one has two scaling variables, or equivalently, two 
moment variables u> and v (see ([)])). The amplitude should now be expanded in terms of double moments 
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with respect to them. There are associated with these double moments new sets of operators that have 
overall derivatives in front (see (^)), which mix among themselves under renormalization via an anomalous 
dimension matrix. The reduced non-forward matrix elements of these operators are thus double distribution 
functions, but they do not seem to have a simple probability interpretation. We will focus our attention 
on the high energy behavior of these distribution functions, solving their renormalization group equations, 
and show that in the high energy limit the gluonic double distribution function actually reduces to the 
conventional (forward) gluon density and the high energy behavior is the same as obtained by a conventional 
double leading logarithmic (DLL) analysis. 

The outline of the paper is as follows. In Section II we define the process and the kinematics of the 
non- forward scattering amplitude T Mt/ . A general tensorial decomposition of into invariant amplitudes 
is given. In Section III we perform a general operator product expansion of T M „, define new moment 
variables and give their relationship to the more conventional Bjorken type scaling variables. In Section 
IV we write down the renormalization group equations, present the equivalence of evolution kernels for the 
double moments, and calculate explicitly the lowest order Wilson Coefficients. In Section V we go to the 
high energy limit to solve the renormalization group equations and make connection with the DLL analysis. 
Section VI gives conclusions and outlook for future work. 



II. DECOMPOSITION OF THE AMPLITUDE 



Consider an unequal mass (non-forward) virtue compton scattering amplitude T M „ as shown in Fig.|j], 
where the incoming photon q' = q — r and the outgoing one q have different invariant masses: 

q l2 = -Ql q 2 = -Ql (1) 



q' q 




FIG. 1. Non-forward Scattering Amplitude 

We restrict our discussion to amplitudes. Cross sections are obtained from the square of the amplitudes. In 
the non-forward case there is not a simple optical-theorem relating the imaginary part of an amplitude to a 
cross section ifTof . Our main interest is in the high energy limit, so similar to || we take all but one of the 
components of the proton momentum to be zero, same for the momentum transfer r. Thus — t = r 2 = 0, 
while r is proportional to the external momentum p and 

r = (p, with C = % ~ g '~ ■ (2) 

2p-q 

We can write the spin averaged amplitude, up to an overall (5-function in momentum space, as 

<?',<?) = i J d 4 ze-^ ( p '|T jM (-|)^(|)b) (3) 

where q — \{q' + q) and —q 2 = \{Q\ + Q 2 ) = Q 2 is now the 'natural' scale of the scattering process. DIS 
corresponds to q = q' while DVCS corresponds to q 2 — 0. We always suppose Q\ < Q\ so that Q > 0. 

Use translation operator one can prove that the conservation of the electro— magnetic current now requires 
the amplitude to satisfy 

^T^ = 0; q' v T^ v = Q. (4) 
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Tu„ is no longer symmetric, rather it is invariant under [i <-> v, q <-> —q'. Together with current conservation 
(Q) , we obtain the general tensorial decomposition of T^ y as 

^=TV = + Tl + (p^P" ~ ^rip^ + i'^Pu) + {^jt) 2 q'^qi^j T 2 , (5) 

where T; = T^g' 2 , q 2 ,p ■ q) are invariant amplitudes analogous to those in DIS and they are even functions 
of p ■ q. We have pulled out from them an explicit factor of \/l — C coming from the external spinors as was 
done in, e.g. , (2). Note we still have only two independent invariant amplitudes even in the presence of three 
invariants, instead of two as in the forward case. 



III. OPERATOR PRODUCT EXPANSION 

In the short distance limit, we can perform an OPE for T Ml/ as a sum of products of local operators and 
their corresponding Wilson Coefficients |ll]] : 

00 J Uj 



J—l n=0 i=l 
00 J u J 



where A^ v and B^ a f} are conserved tensor structure operators and all indices except /i, v are symmetrized 
as usual. B^ap is the same as in the forward case, i.e. , 

B^vap = g^agvp^ + g^vdadp — g^udvdp — g v pd^d a , (7) 

and corresponds to the tensor structure multiplying T2 in (Q) while we have not found the explicit form of 
A^v which would generate the tensor structure corresponding to Ti. But, as in the forward case, once we 
know T 2 , we can obtain T x by a Callen-Gross relationship, at least in leading logarithmic level (see Jl6|)). 

Because now the momentum flowing into the local vertex is r instead of zero, we have to include in the 
expansion new sets of operators that have overall derivatives. In leading twist (the meaning of which in this 
case will be clear later) the operators for QCD are 

tylffl = d^...d^ g(0) 7M „ +1 i V» n+2 ...i g(0) (8a) 

= d^.-d^ F^(0)i V^ +1 ...» V„j F^(0). (8b) 

After taking the matrix elements between the asymmetric external states, and after taking a Fourier trans- 
form, the external derivatives would eventually be turned into factors of r M while the internal derivatives 
into either r M or (p + p')^ = (2p— r) M . We thus define the two moment variables in the non-forward case as 

(2p - r) ■ q r-q 

The forward case would be the limit v — while DVCS corresponds to v — 1. 

In a QCD-parton picture, the diagrams contributing to T M „ are the so-called 'hand-bag' diagrams as shown 
in Fig.|^ If we parametrize the momentum of the scattered parton as k — xp + yr (see, e.g. |2j), where x 
and y are two Bjorken type scaling variables conjugate to u> = 2p- q/Q\ and v = 2r ■ q/Q\, respectively, our 
moment variables are related to these quantities via 

w = %(xu> + yi>- \i>) ; v = (10) 

Q 2 Q 2 
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q' 



q 



q' • y '\ 



xp+yr 




xp+yr 




p p \ / p p 

FIG. 2. Handbag Diagrams in QCD-Parton Picture 

There will be mixing between operators with same J but different n labels. To be more specific, evolution 
in Q will lead to mixing from 0^ J ' n ^ to 0^ tyJ ' n ) with < n < n! < J. But we do have the freedom to 
choose, for simplification, at the factorization scale no, that all internal derivatives give factors of (2p — r). 
Thus we can write 



and the reduced matrix elements will depend only on J — n: 

(p'||6»W")||p) = (p'||6«(- 7 -"'°)|b) . 



(ii) 



(12) 



Analogous to the forward case, taking the logarithmic derivative of q 2 of the Fourier transformed Wilson 
Coefficients E wc define the Wilson Coefficients in momentum space E as 



Q 

where we have pulled out an explicit factor of =| to make the form of E simple (see J3C|)). We have 

Q 



(13) 



T^=e yj<p'i|6w«")||?>) 

J,n,z 

l 



Mo 



p-q, , / ^ , (P-i-,2 , \ p(.) / w \ 

. P^Pv - —{P^qu + q^Pv) + (—r) %q. v ^j_2,„(a s , — ) 
P'Q V 9 9 / Mo 



(14) 



F' is a combination of £ and a similarly defined F. Note we have explicitly indicated their dependence 
on the factorization scale. Non-leading twist terms in this case are suppressed by a power of =2-. Thus we 

obtain the invariant amplitudes as follows: 



Q 



J,n,z 



p-<? 

M 2 



J,n,i 



VWT 2 = - £ <p'||O«(^|b) (p0) c^V"F« 2! Ja s ,-^) . 

Mo 



(15a) 
(15b) 



IV. RENORMALIZATION GROUP EQUATION AND ANOMALOUS DIMENSIONS 

A. Renormalization Group Equation 

We concentrate on T 2 . Ti can be obtained, as in the case of a forward scattering, from the Callen-Gross 
relationship 
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Ti = ^T 2 (16) 
T2 can be regarded as a double distribution function and is the sum of 'double' moments 

T 2 = ^w J - r V l T^' n) , (17) 

with the double moments T^'^ as 

^IXflL) =-E(P / |P (i)(J '" ) lb)(,o)^ 2 ,n(^>f ) (18) 

We will eventually analytically continue in J, but leave n as discrete. The inverse of expansion ( |l7j ) is 

J c U) J c V UJ V 

= 2i r^/^ ( i)-(W,^) (19) 

Ji u J c v u v 

where c is any contour in the io(y) plane that encloses the origin and W2 = I111T2. Note in obtaining the 
second equality we have used the fact that T 2 is even inw. T 2 now as a double distribution function is given 
by 

T 2 (Q 2 , uj, v)= f^~J2 eV-^ l °z" v n T 2 J ' n) (20) 

with the contour c lying parallel to the imaginary axis in the J — n plane and to the right of all singularities 
in that plane. 

Under renormalization, these operators scale according to a renormalization group equation 

M 2jLo««n) = J- ^rLO^^ (21) 



n —n 1 



where T is the anomalous dimension matrix which is upper triangular and acts in the product space J ®uj. 
The Wilson Coefficients then obey a similar renormalization group equation 

f^&^M, ? } = £ X^S f ™ ? »- (22) 



We can write the solution as 



frM-^EE^S^ 1 ) M# n (a(Q),£), (23) 



where M is a path ordered exponential of the anomalous dimension matrix, formally given as 



Mtf n (a(Q),Q)= (veM-J 2 f(a s (A 2 ))^)j (24) 



Thus the double moments are given by 



^fCj'-EI) 1 ^^^ m£>(Q),£) d/lldww^iw^. (25) 
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B. Anomalous Dimensions 



In light-cone gauge, the lowest order anomalous dimensions of the quark and gluon operators are generated 
by the usual triangle diagrams (note: the graphs where we have more than two quark/gluon lines meet at 
the vertex are zero because of the gauge choice) as shown in Fig. || 




-qq Igq Iqg Igg 

FIG. 3. Triangle Graphs for the Anomalous Dimensions 



Compute the graphs using light-cone variables with the definition 

w=— , v=— , (26) 
fci+ k 1+ 

and the conventions (cf. |fl2|| , also see Appendix A) 

0[ J ' n ^ + \k, k-r)= 7+(2fc - rY^rl (27a) 
(J,n)*f3( kj k _ r j = 2ga/3 n .kn-(k-r)(2k- ryl^r'l (27b) 

where n 2 = and n ■ A — A+ for any four vector A, yields (see Appendix A for details) 
/7r Jo 



(2w - z/) J ^ S(v-u) -u(l + 



V 1 — U! 

Q(cj - v) J- 1 + C " 2 ~ + * (1 - w) (-2I - log(l - 1/) + ~) ) ( 2N0 

1—vl — w 2 ' 



rM=^<7 F7 , fcfrV / ^(2^-z/)- 7 ^ 1 - 2 f e(^-cj)i(w 2 -2w)-e(cj-^)^(i+(i-w) 2 -^)^ , (28b) 

27T Jo V v 1 ~ U J 

Ki = ~2g ij k(,v n [ diu{2j~^ l - 1 (^^~Lu)-uj(2ij-l-u)+e(io~^)^{uj 2 + {l~Lo) 2 -^ , (28c) 

27T 2 Jq V f 1-1/ J 

T r, a^A 2glJ k j n f 1 ^ {2uj _„ )M (q (1 (4^^ +M + ( " 3+ " 2 ' 2 " 2 " -3(c 2 +c))) 

MJ 27T Jo V ^ f 1 — 

l — i/ 1 — cj 1 — w 2Ga / 



We have included virtual corrections coming from the self- energy graphs that will cancel the colinear singu- 
larities at end point in the u> integration (Jo terms, see (A9) for the definition of Jo) and give correct constant 
terms in the anomalous dimensions. is the usual step function and &o — H — f n / is the leading coefficient 
of the QCD /3-function. 

Some comments are in order. It is straightforward to show that in the forward limit where v = 0, ( p8| ) 
reduces to the conventional Altarelli-Parisi splitting functions However, because of the no n- forwardness, 
there is no simple probability interpretation for these evolution kernels as splitting functions. 

The reason we have two terms for each triangle graph is that for different fc+ integration regions pinching 
of the fc_ pole is different (see Appendix A). Thus there is no straightforward optical-theorem type dispersion 
relationship between the cross section and the imaginary part of the amplitude. Nonetheless, we can still 
do analytic continuation in J and relate the matrix elements to 'double parton distributions' flOl which now 
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may not have a direct probability interpretation. The second terms of our T differ from those of [|6| only 
because the factor n ■ kn ■ (k — r) in our (p7|b) is not included in the definition of their gluon vertices. 

After we perform the to ( i.e. k + ) integral, we can in principle obtain the corresponding anomalous dimen- 
sion matrix in moment space. At first sight, the (1 — v)^ 1 factor of the second terms may generate a series of 
infinite sum over powers of which will spoil the locality of the vertex thus invalidate the operator product 
expansion. Detailed calculations show that because the lower bound of integration is now v instead of zero, 
there will also always be at least one power of (1 — v) coming out from the integral. Furthermore, the highest 
powers of v cancel completely between the first two terms of each T, leaving the highest surviving v power 
as i/ J_1 for the quark sector and v for the gluonic sector, as must be the case to make the OPE valid. All 
the Io and log(l — v) terms also cancel completely between real and virtual graphs. 

Because of the mixing among different n moments, as well as between quark and gluon sectors, it is very 
difficult to read off the anomalous dimensions of the mixing between two operators (with same J) of definite 
n moments. However, the form of the anomalous dimension matrix simplifies when the high energy limit is 
taken, and we will be able to make connection between this general formalism and the conventional leading 
logarithmic approximation (LLA) analysis results. 

C. Lowest Order Wilson Coefficients 

The lowest order Wilson Coefficients can be calculated in perturbation theory from tree level diagrams 
like the upper part of Fig. |^. The result is 

00/ / 
t(0) a ' 2 fi 7V^ 2i+i I P'Q, , 1^1" \ 

1=0 \Q 1 
+ p~q ~~ + + ^^^^"Y) ■ ( 29 ) 

Comparing to ([l4]) we obtain 

Ej q n°l(a s {Q), 1) = Ue% J odd, J > 3 (30a) 
£^ (0) (a s (Q),l) =0, J, n otherwise. (30b) 

This means at leading order there is no dependence on the non-forwardness in the Wilson Coefficients. 

V. HIGH ENERGY LIMIT 

To exactly solve the renormalization group equations, one needs to diagonalize the anomalous dimension 
matrix in both flavor (i) and moment (n) indices. In the high energy limit, however, the situation simplifies. 
After analytical continuation in J, at high energy, the dominant contributions to T 2 come from the leading 
(right most) poles in the J-plane, while in flavor space t he g luon anomalous dimension dominates. In the 
terms of the expansion of the path ordered exponential (pi|), all the factors of the product of anomalous 
dimensions are V" except the first one, which should be Tis due to the quark loop, i.e. we have to force 
one gluon-quark transition at the end of the evolution. The two flavor sums in (P5|) then collapse to i' = q 
and i = g. We evaluate M between a quark and a gluon state to obtain, again formally, with 7" now the 
anomalous dimension matrix in the moment space, 

M«„(a(Q),^)= (- [ Q ^V(a s (A 2 ))") (Pexp(- ^ 7 ^(a s (A' 2 )))") . (31) 

The leading pole terms of the relevant anomalous dimensions can be parametrized as (suppressing the J 
label) 
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1 4- c qa 

7^(A 2 ) = -^a s (A 2 ) 7» with 7" = ; (32a) 

7^(A 2 ) = -^« S (A 2 ) 7» with 7 f„, = • (32b) 

while < n < n' < J for 7 ran » and c l n l n vanishes when n' — n. Because of the upper triangular structure 
and the leading pole positions of these 7 matrices, we have, since m < I < n, 

^ ^fml Tin = 'Ymn 'Ynn — Tmn a n- (33) 
I 

This means in a product of 7 matrices we can, in leading logarithmic approximation (LLA), take only 
diagonal entries in all the factors except the first one. Thus the second factor on right hand side of ( pl| ) 
becomes 

I + \{e A ^< - 1) 7 f„ 9 = \e A{>?) < 7 f„ 9 + non- leading-terms. (34) 
where A(\ 2 ) = ^ i og Thus, © reduces to 

Q n f f Q2 dX 



M™, n (a(Q), f Q ) = gJ 2 ^« S (A 2 ) >- y«J n . (35) 



The n' summation in the expression for the double moments ( |25| ) also collapses in LLA to the diagonal 
element of j qg , and we have the final expression as 

KvT^WSL , = -BSL.IO.. 5) wia«(«iw u y-_ r d 4«.(*) (36) 



The invariant amplitude T2, after analytically continuation in J becomes 



M 2 



,/WT 2 (qW)= / £ £V /" ^«,(A 2 )^ 

■41,„K 5 £) (p'||0 (9)W " ) lb> w ^exp((J-n)log W + 7 ^) . (37) 



A saddle point approximation in high energy limit for the J integartion leads to 

p ' q P\ 7t n f I 71 -log" dX2 „ f^,/ x 2^i/4 1 

-VWT2 = r t , 3/2 e /. -tt-Q! s (A )A(A ) 



M 2V s ' 2^Vlog 3 / 2 ^ y^g A a"»vw v > Ja^x 

1 * ' log LO 



■eV^ 1 ^ v n E { j% n (a s , — )(p'||6 (9)Wn) |b) ((llo) . (38) 



Because the saddle point fixes the value J — n = J s — n = 1 + y , which fixes the number of internal 

derivatives inside the gluon operator, recall ( |l2| ) it is clear that the reduced matrix element in (§8|) is the 
conventional (forward) gluon distribution function in moment space with a shifted moment label J s — n. 
Despite the formal summation over n, which is essentially the only difference introduced in this limit by 
the non- forwardness, only one non-perturbative input is needed, which is the gluon density. This means 
that although in general there will be complicated dependence on v in the double distributions, (in high 
energy and hard scattering limit) it is nontheless perturbative. At leading order the n summation collapses 
when we take the lowest order Wilson Coefficients given in (30) with the flavor averaged quark charge 
nf < e 2 . >= J^f e q ■ The final expression of the double distribution T2 is 
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P'9 fT^rv in 2 \ m f <e l > { 1 I 

f^ TjW ^ )= "v^r log5 ZW) ) 

■{p'\\6^ J ^\\p) M e l °e« e V^ A( ° 2)logw . (39) 

We can see clearly that the leading high energy behavior is exactly the same as that given by a forward 
direction leading logarithmic analysis. 



VI. CONCLUSION AND OUTLOOK 



We have formulated a general operator product expansion of a non-forward and unequal mass virtual 
compton scattering amplitude. We found that, because of the non zero momentum transfer, the expansion 
now should be done in double moments with respect to the moment variables defined in (|9j). The double 
moments can be parametrized as products of Wilson Coefficients, which can be computed perturbatively, 
and non-forward matrix elements of new sets of quark and gluon operators, namely, double distribution 
functions. These operators mix among themselves under renormalization and they obey a renormalization 
group equation. We calculated the equivalence to evolution kernels of double distribution from which the 
anomalous dimension matrix of the operators can be extracted. 

In high energy limit, we found the leading contributions of the anomalous dimension and used them to 
solve the resulting renormalization group equation. We recovered the conventional DLL analysis results 
and made the connection that in fact the double distributions are proportional to the conventional forward 
gluon density in this limit. This justifies previous analyses using DLL and forward gluon density on non- 
forward processes like the exclusive vector meson production (e.g. H). Furthermore, we have the formalism 
to proceed in principle beyond DLL to next to leading order. 

Our analysis is quite general, we can go from DIS (forward case, corresponding to v — 0) to DVCS 
(corresponding to v = 1), while diffractive vector meson production {y close to 1) can also be studied by 
replacing the second quark-photon vertex with effectively a wavefunction of the vector meson. 

There are still more details of the analysis that need to be filled in. The explicit form of the anomalous 
dimension matrix is not yet written down, the dispersion relationship between the invariant amplitudes TVs 
and the double structure functions need to be clarified (for example, the exact physical meaning of the fact 
that different fc_ poles are being picked out at different fc+ integration regions). Some work towards these 
ends is still in progress. 
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APPENDIX A: 



In this appendix we will show in some detail the computation that leads to (|28|), in particular, the quark- 
quark anomalous dimension. 

Usingstandard Feynman rules and light-cone (LC) gauge we can write down the value of the first diagram 
in Fig. H as 

where the quark vertex O q is given in ( p7[ ) and the LC gluon projector (numerator of the LC gluon propagator) 
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D"P(k) = g"? - + , fc " n/3 (A2) 
n ■ k 

with n 2 = and n ■ A = A+ for any four vector A. Writing the integral of the loop momentum in terms of 
light-cone variables we get 

99 (2tt) 4 7 _ + (fc 2 + ie)((fc - r) 2 + «e)((fci - fc) 2 + le) v y 

where 

A(fc,r) =7 a (/!- /) 7 + MpD^ih -k). (A4) 

We perform the fe_ integration first, in terms of the poles of the integrand. However, the positions of poles 
in fc_ arc different when k + is in different regions. More specifically, when k + < and fc+ > fci+, all the 
poles in /c_ lie at the same side of the real fc_ axis and we can distort the contour so the integration over 
fc_ gives zero. On the other hand, when < k + < r + , we can complete the contour in the lower-half-plane 

and pick up the pole at fc_ = -M ie = k^\ while when r + < k + < k\ + we complete the contour in the 

k 2 (2) 

upper-half-plane to pick up the pole at fc_ = 2{k~-k 1+ ) + * e = ^- • After evaluating the residues we have 



_ _ t ^L ~ 2ni f r+ dk^A^(2k - r ) M r n 



2 "' ' 1+ dk+ih- k) + A ( - 2) (2k^r) J ^ ri - 1 rl J , (A5) 



-2fci+(fci - r) 



+ 



+ 



where A^ and A^ 2 ) are A(fc,r) evaluated at the poles k^ and k^\ respectively. We are interested in the 
logarithmic divergence in the graphs so in computing them we keep only the leading powers of k which is 
quadratic. We obtain 

A™ = -2k 2 1+ (l + — — ) (A6a) 

1 — UJ 

A(2)= -(Sp (1+ " 2 - i/(1+ " )) (A6b) 

with uj and v as defined in (ffit). Thus we have 

C" - a fo, J § 7 + *frV (i jf *- (2- - » (i + ^ 

' ^ (2„ - ,,)•<-' 1 + " a -"' 1+ "» N ) . (A7) 



\ — vJ v 1 — UJ 

The u; integration is divergent at u> = 1. This divergence is cancelled by the self-energy graphs, the value of 
which in LC gauge can be readily taken from || as 

Z F ( X ) =1 + 2^1 (-27 - 21og H + (A8) 

Z7T e z 



where Iq is the colinear divergence 



1 dz f 1 dz 



Io= ~= , (A9) 

Jo z Jo 1 - z 



and Zp depends on the longitudinal momentum fraction x. So in our case the self-energy contribution from 
the k\ — r line should have x — n ■ (k± — r)/n • ki = 1 — v while the k\ line has simply x = 1. Adding 
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^(Zp(l) + Zp(l — v)) to (A7) and identifying the logarithmic divergence J* ip- = ~ we will obtain the 
expression of T qq as in (p8|) . Note there are no other diagrams contributing in LC gauge. 

The computation of the gluonic sector follows similarly. The one point one needs to pay attention to is 
the form of the gluon vertex in (B7| ) . The full tensorial structure of the LC gluon vertex in the non- forward 
case, generalized from that in [jljjjis 

K*a(^' k — r) = gapn-kn- [k — r) — n- [k — r)n a kp — n-k(k — r) a np + k- [k — r)n a np. (A10) 

This vertex is to be contracted with gluon lines D aa (k) and (k — r). It is straightforward to verify that 

D aa '(k)V^ 0l (k,k-r)D^'^{k-r) =V^(k,k-r) (Alia) 
D aa ' {k) g a > p> n-kn-{k ~ r) D (i ' (i {k - r) = V° p {k,k-r), (Allb) 

thus we can substitute the full gluon vertex by its first term because of the LC projectors of the two gluon 
lines connected to the vertex and use (E7 



X. Ji, hep-ph/9609381, to appear in Phys. Rev. D . 
A.V. Radyushkin, Phys. Lett. B380, 417(1996) 
S.J. Brodsky, et al. Phys. Rev. D50, 3134(1994). 

A.V. Radyushkin, Phys. Lett. B385, 33 3( 1996) 

J.C. Collins, L. Frankfurt, M. Strikman, }iep-ph/9611433| 

L. Frankfurt, A. Freund, V. Guz ey, M. Strikman |hep-ph/970344l 
X. Ji, W. Melnitchouk, X. Song. |hep-ph/970237g|. 



A.V. Radyushkin, ^iep-ph/9704207 . 

G. Curci, W. Furmanski, R. Petronzio, Nucl. Phys. B175, 27 (1980). 

N. Christ, B. Hasslacher, A.H. Mueller, Phys. Rev. D6, 3543(1972). 

J. C. Collins, Renormalization, (Cambridge University Press, London, 1984). 

A.H. Mueller, Phys. Rev. D. 18, 3705 (1978). 

G. Altarelli and G. Parisi, Nucl. Phys. B126, 298(1977) 



11 



